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Spin dynamics as well as static properties of the one-dimensional J-J' model (S = 1/2, J > 
and < a — J' / J < 0.5) are studied by the exact diagonalization and the recursion method of 
finite systems up to 26 sites. Especially, the dynamical structure factor S(q, u>) is investigated 
carefully for various values of a. As a increases beyond the gapless-gapful critical value a c = 
0.2411, there appear features definitely different from the Heisenberg model but the same with 
the Majumdar-Ghosh model. Some of these features depend only on the value of a and not on 5: 
a parameter introduced for the coupling alternation. By comparing these results with a recent 
inelastic neutron scattering spectrum of an inorganic spin-Peierls compound CuGeOs [M. Arai 
et al. : Phys. Rev. Lett. 77 (1996) 3649], it is found that the frustration by J' in CuGe0 3 is 
unexpectedly strong (a = 0.4-0.45), and at least a must be larger than a c to some extent. The 
value of J is evaluated at ~ 180K consistent with other estimations. The coupling alternation is 
extremely small. This large frustration is a primary origin of the various anomalous properties 
CuGeOs possesses. For comparison we refer also to a'-NaV^Os. 
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§1. Introduction 

In connection with the quasi spin gap found in the 
high-T c superconductors, materials with a spin gap have 
been intensively studied both by experiment and theory. 
For a systematic understanding of the spin gap, it is in- 
dispensable to investigate thoroughly the properties of 
basic models. Among such models, the 5 = 1/2 one- 
dimensional (ID) antiferromagnetic (J > 0) Heisenberg 
model with frustration by the next-nearest-neighbor ex- 
change, which we call the J-J' model,EP 



H = j'£s r s j+l + j"£s j -s. 



3+2, 



(1.1) 



has some unique points. It is well-known that this model 
changes its character in the low-energy excitation from a 
gapless feature of the Heisenberg model (a = J' j J = 0) 
to a finite-gap excitation like the Majumdar-Ghosh (MG) 
case (a = 0.5) at a c = 0.241 iff At the MG point the de- 
generate ground states are exactly known as .products of 
the nearest- neighbor singlet wave f unctions. a cP So far, 
some static and thermodynamic properties have been 
elucidated by numerical studies™ and low-energy proper- 
ties by the conformal field theory.aEP However, system- 
atic studies of spin dynamics have been lacking™ 

Another important aspect of this model is a close re- 
lationship with the first inorganic spin-Peierls (SP) com- 
pound CuGe0 3 .Q) The origin of the spin gap in CuGeOs 
is now established as the lattice dimerization,Er EP al- 
though it is very small. Some properties below the SP 
transition temperature Tgp ^14K are common with the 
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organic SP compounds and can be. broadly understood 
by the theory of Cross and Fisher S3 On the other hand, 
there are many unique aspects contradictory to this the- 
ory. For -instance, phonon softening has not yet been 
observed, 113 which has been necessarily discovered for 
the organic SP compounds as a driving force of the SP 
transition. Moreover, the thermodynamic properties ob- 
served above Tsp do not coincide with the results of the 
simple Heisenberg model .cP 

Later, it was found that the latter contradiction can 
be reconciled by introducing frjistration due to the next- 
nearest-neighbor exchange. ^ZLaES Although many fol- 
lowing theoretical studiesEj'ESEZP have unanimously con- 
firmed the importance of this frustration to explain the 
experiments, the ratio of frustration a remains contro- 
versial among them. Originally, Riera and DobryEj con- 
cluded a — 0.36 by fitting_^the spin susceptibility x(T)- 
Meanwhile, Castilla et ality excluded the possibility for 
a > a c by a conjecture based on the theory of Cross and 
Fisher and the dependence on temperature of the gap, 
and concluded a — 0.24, which is a marginal value for 
the gapless condition. At any rate, the main point of the 
controversy is how the gap by J' for a > a c is compatible 
with the SP transition. 

In this paper, first we investigate spin dynamics of the 
ID J-J' model without coupling alternation carefully for 
various values of a (0 < a < 0.5) by the exact diagonal- 
ization and the recursion method at zero temperature. EiP 
As a result, the dynamical structure factor S(q,u>) has 
qualitatively different characters, according as a < a c or 
a > a c ; for a ~ a c it resembles that of the Haldane- 
Shastry (HS) model. 

Second, we compare these results with a complete 
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spectrum of the inelastic neutron scatteringB to search 
proper parameters extensively in the a-5 space, by 
adding the effect of coupling alternation 5 to the J-J' 
model eq. (1.1). Note that our fit is carried out at low 
temperatures 7sp) in contrast with the previous 

ones by x(^)Ji3't-J which were done for T > 55K. Con- 
sequently, we have found that a thus obtained is unex- 
pectedly large (a = 0.4-0.45), and that at least a must be 
somewhat larger than a c . This conclusion is drawn from 
the fact that the experimental spectrum shows some ev- 
ident features characteristic of the region a > a c , which 
do not appear only by adding 6 to eq. (1.1) with a < a c . 

Once a is specified, the value of J is determined by 
neutron dataEf as ~ 180K, ffihM 1 i s consistent with 
those obtained by other meansiiirlEP On the other hand, 
5 is found to be very small (0.001-0.005). We are con- 
vinced that the main origin of most unique properties in 
CuGeOs is nothing but this large frustration, which is 
not in competition but in concert with the lattice dimer- 
ization in making a gap. In this connection, we refer 
also to the second inorganic SP compound a'-NaV^Os, 
in which a seems small. 

This paper is organized as follows: In §2 we mention 
the model and the method. Section 3 concentrates on 
the results for the J-J' model without coupling alterna- 
tion (J = 0). In §3.1 the Heisenberg case is taken up 
as a reference. In §3.2 the regime of a ~ a c is consid- 
ered in relation to the HS model. In §3.3 characteristics 
of the gap regime are explained. Section 4 is assigned 
to the comparison with the experiments of CuGeOs. In 
§4.1 we compare S(q, u) with a complete spectrum of the 
neutron scattering to determine the value of a. In §4.2 
the values of J and 5 are considered. In §4.3 based on 
the parameters thus determined, we discuss a couple of 
other properties observed in CuGeOs and a'-NaV~205. 
We summarize in §5. In Appendix A we describe the 
single-mode approximation for the J-J' model. In Ap- 
pendix B we give a proof of a pure single mode at q — it/2 
for the exactly soluble case: 2a + 5 = 1. 

§2. Model and Method 

In this section we summarize our model and method 
briefly. In this paper, we consider the ID antiferro- 
magnetic Heisenberg Hamiltonian with the next-pearest- 
neighbor exchange and the coupling alternations 

n = J E{ t 1 - S J ■ + aS i ■ S i+*}> 

3 

where we assume a, 5 > and Sj+N = Sj with N being 
the number of sites and even. 

For a — 5 — eq. (2.1) is reduced to the Heisen- 
berg model, for which various ground-state and ther- 
modynamic properties are exactly known by the Bcthc 
Ansatz; quite recently the structure of spin excitation_at 
T = has been clarified by using quantum group, etc.c$ 
The ground-state phase diagram of eq. (2.1) in-the a-5 
plane was given by the conformal field theory pi? the 
regime of gapless excitation is limited on the segment 
< a < a c and 5 — 0. At a c marginally irrelevant op- 
erators vanish (conformal invariant) , so that there exists 



no logarithmic correction. The critical value was deter- 
mined as ~ 0.24110) by using the fact that the lowest 
singlet and triplet excited levels cross at this point in 
the thermodynamic limit. For a c < a < 0.5 and 5 = 
an excitation gap opens as in Fig. 11(b), and the ground 
states are singlet and twofold degenerate in the thermo- 
dynamic limit. But this degeneracy is lifted for finite N 
and a ^ 0.5. 

At a. pj= 0.5 the degenerate exact ground states are 
knownEftf as: ^± = P±(ipi ± V'2), where 

^i-[l,2][3,4]--.[A-l,A], 

^ 2 = [2,3][4,5]---[A,l], (2.2) 

and 



v± = 1 



/ ^2 N / 2 + 1 ± (-l)^/ 2 



Here, the square bracket indicates the singlet pair spin 
function: 

[i,j}=a(i)(3(j)-PdMj)- (2-3) 

*+ (v?_) belongs to Q — (tt)& The spin excitations 
for the MG point wepe studied by introducing an varia- 
tional wave function,© in which an isolated doublet (or 
spinon) pair plays an important role. 

For other values of a in the J-J' model eq. (1.1), the 
ground-state and thermodynamic properties are stud- 
ied by the exact diagonalization and a Quantum Monte 
Carlo method.H' 

For a finite value of 5, a finite gap always exists. _The 
exact ground state is known on the line 2a + 5 — 1 £P In 
this case, ipi, in which singlet pairs sit on the bonds with 
the larger exchange J(l + 5), becomes the unique ground 
state, and ip2 (therefore ^f±) is no longer an eigenstate. 
For 2a + 5 > 1 the phase is "spiral" . In this paper, we 
concentrate on the parameter range 2a+5 < 1 (a, 5 > 0), 
and leave the spiral phase for future studies. 

For this model, our main concern is the spin dynamical 
structure factor at zero temperature: 



S(q,uj) = ^\(^ n \S^ )\ 2 5 (u; - (E n 



Eo)), (2.4) 



where fi — 1, $„ (E n ) is the n-th eigenfunction (eigen- 
value) of the system, and the ground state is indicated by 
n = 0. n runs over all the eigenstates. S(q,u>) is a quan- 
tity proportional to the spin contribution in the intensity 
of an inelastic neutron scattering spectrum. Since here 
the ground state is always singlet for finite systems, only 
triplet states contribute to S(q,u>) due to the selection 
rule. For finite systems, S(q,u>) can be estimated by the 
recursion (or continued-fraction) method.tSP To this end, 
we expand the corresponding Green's function G(z) as 

G{z) = (* |^ 9 _2_^|vl/ ) = S(q)C(z), (2.5) 
where S(q) = (* |SI „^|*Q> and 



C{z) = 



z — ao 



(2.6) 



z — a\ — 



Z — 0,2 
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For finite systems, C(z) is comprised of a finite number of 
S functions. Given the exact ground-state wave function, 
which we usually obtain by the exact diagonalization, 
we can determine the coefficients a^s and 6,'s through 
a Lanczos-like procedure. Using this expansion, we can 
rewrite eq. (2.4) as a sum of residues, 



S(q, <jS) = lim ImG(uj 

= S(g)Re^Res 



if] ~ Ea) 
C{LJ + ir)-E Q ),E^\, (2.7) 



where t runs over all the poles of C . 

Here, we write down some technical points. Actually, 
we calculate the contribution to S(q, u>) by integrating C 
for each pole with finite r), typically 1-5 x 10~ 5 J. Since 
the sum rule 



S(q) = I dui S(q,ui) 
Jo 



(2.8) 



holds, the total contribution of residues for a fixed value 
of q is normalized to unity irrespective of the value of 
7/. We have checked that the total contribution in our 
calculation is always ^ 0.9999. Next, to raise preci- 
sion, it is important to start with a ground-state vector 
as precise as possible. Otherwise, numerical errors are 
rapidly mingled in the coefficients. In fact, a^s and frj's 
accurate enough to be used seem to be limited to i ^ 20- 
30 at most, so that we cut off higher coefficients than 
i = 30 mostly in this study. This hardly affects poles 
and residues in low energy and with strong intensity, but 
causes inaccuracy to some extent of the ones in high en- 
ergy and with weak intensity. In this context, however, 
the results are justified for most cases by confirming the 
regularity of the dependence on system size. 

Finally, we mention the finite-size scaling. For the J- 
J' model, there is no simple scaling function which fits 
finite-size data correctly.cP However, by searching for a 
pertinent function case by case, one can fulfill an extrap- 
olation reliable enough to discuss the physics. We use 
the usual polynomial fit of order.-smaller than 10th, in 
addition to the formula for 5 ^ 0x2$ 

N ' 



r(jv) = r(oo) + 



7 

N rn 



■ exp 



N 



(2.9) 



where 7 and No are fitting parameters, and m — 2 (1) is 
used if r is energy (a quantity of energy difference like 
the gap A). 

§3. Spin Dynamics of the J-J' Model 

In this section, we focus on the J-J' model eq. (1.1) 
without coupling alternation. By contrasting these re- 
sults with the finite-5 cases, the effect of a is definitely 
discriminated from that of 5. 

3.1 Excitation spectrum of the Heisenberg model 

Although the spin excitation of the ID Heisenberg 
model has been studied for a long period, it was quite 
recently that the accurate behavior of the spectrum was 
clarified. In contradiction to the spin-wave picture (spin 
1), the excitation spectrum of the Heisenberg model con- 
sists of even number of free-spinon (spin 1/2) excita- 



tions.^ As the lowest-order contribution, two-spinon ex- 
citations form a dominant continuum, the boundaries of 
which coincide with the known curves!^ 



, , nJ . 



and uj u (q) 



7r J sin — . 



(3.1) 



We call this continuum as the two-spinon continuum 
(TSC) in this paper. Some years ago, an approximate 
formula of S(q, to) for the TSC was proposed by Miiller 
et al. from the consideration of the exactly soluble-Jf y 
model, finite-size results and various sum rules, as:eZP 

S{q, to) = f = Q(u - ^)6K - uj). (3.2) 

Here, A is a constant of order 1, and O is the step 
function. Actually, the intensity of this spectrum was 
roughly reproduced by the neutron scattering experi- 
ments of KCuF3^ a good ID antiferromagnet. On the 
other hand, they simultaneously pointed out that it de- 
viates slightly from the exact one mainly due to higher- 
energy excitations outside the TSC. 

Recently, the two-spinon and higher-order-spinon 
(HOS) contributions to S(q,co) have been exactly calcu- 
lated by using formulae of quantum groupjHj and com- 
pared with the pcact diagonalization as well as the above 
approximation]^ According to these studies, although 
the formula eq. (3.2) approximately represents the ex- 
act two-spinon spectrum, eq. (3.2) misses the behavior 
in the vicinity of both boundaries. Especially, it overes- 
timates near the upper boundary, under which there has 
to be a square root singularity in S(q,ui), instead of a 
jump given by eq. (3.2). Meanwhile, the two-spinon con- 
tribution to S(q,cu) occupies 72.89% of the total weight; 
the HOS processes have not a little effect on S(q,u>) as 
expected before by Miiller et al. 

Here, we summarize the characteristics of the spec- 
trum of the Heisenberg model as a reference for later 
discussions, based on our calculations by the recursion 
method. We note that there have been a number of 
studies in this line for the Heisenberg model. In Fig. 1 
S(q,ui) of the Heisenberg model is shown for N = 26. 

Poles by the two-spinon processes have dominant in- 
tensity and are situated on a series of sinusoidal curves 
in the TSCEf The characteristics of the TSC are [a] 
the intensity becomes strong as q approaches it due to 
the logarithmic divergence of S(q) at q = ir. [b] At ev- 
ery fixed value of g, S(q, u>) is monotonically decreasing 
function of u>. 

On the other hand, the other poles are of the HOS con- 
tribution, which continuously spread inside (q/ir ^ 0.4) 
as well as outside of the ISC. The existence of this contri- 
bution was suggestedEJ'EZP in relation to the single- mode 
approximation as discussed in Appendix A. The inten- 
sity of these poles is considerably weak in comparison 
with the two-spinon contribution, but severely depen- 
dent on N and abruptly increases, as seen in Fig. 3(a). 
The lower boundary of this HOS continuum is indicated 
by dashed lines; it seems the lowest mode switches a 
couple of times, and this boundary lowers with increas- 
ing N. Incidentally, this curve looks similar to that 
of the lowest quintet (S — 2) excited levels (solid dia- 
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mond forJV = 14), which are described by four-spinon 
processesBiP The four-spinon continuums of the triplet 
states and of the quintet states may coincide in the ther- 
modynamic limit, just as the TSC's of the triplet states 
and of the singlet ones do.EP At any rate, another im- 
portant feature in the spectrum of the Heisenberg model 
is [c] there is non-two-spinon contribution, although its 
intensity is weak. 

In the remainder of this subsection, we consider the de- 
pendence on system size of the poles and their residues to 
discuss whether a pole belongs to a continuum or forms 
an isolated branch. Such a finite-size analysis was suc- 
cessfully used for the detection, of an isolated branch in 
the S = 1 Heisenberg model,E3' although we cannot al- 
ways draw a definite conclusion from this analysis. 

Following this work, we plot residues of the poles be- 
longing to the lowest two two-spinon branches in Fig. 
2(a). Since both branches belong to a continuum, the 
value of each residue ought to vanish in the thermody- 
namic limit. As fox-.the branch 1', the weight decreases 
with increasing NS3 However, the weight of the branch 
2' is almost independent of N for large q and rather 
increases with increasing N for medium q within these 
system sizes. In fact, we have checked various cases and 
found that the dependence on N of residues is sometimes 
not monotonic. 

On the other hand, Fig. 2(b) shows the dependence on 
N of the positions of poles for the lowest four two-spinon 
branches. Although the branch 1' is scarcely dependent 
on N, other branches have considerable dependence on 
N. In the inset of Fig. 2(b) the pole positions at q = ir 
are plotted versus 1/7V 2 . All the branches monotonically 
decrease as the system becomes large, and seem to gather 
at oj = 0. Note that S(q, lo) diverges as l/oo at q = tt. 

As seen in this case, a pole which belongs to a con- 
tinuum tends to have appreciable dependence on system 
size at least either of its position or of its residue. On the 
other hand, if a pole has little dependence on N both of 
its position and of its residue, and unless adjacent poles 
converges to its position, that pole probably forms an iso- 
lated branch. However, as will be discussed in §4.3, this 
will not be a necessary condition for an isolated branch. 
Thus, one should check carefully the dependence on N 
of the positions and residues in each case. 

3. 2 Excitation spectrum near the critical point a c 

As mentioned in §2, at a c = 0.2411 the system eq. 
(1.1) is conformal invariant. This fact reminds U£Lof an- 
other conformal invariant model, the HS modeller 

n = %W^W Si ' Sh (3 ' 3) 

where d(r) = (N/ir) sin(7rr/A). This model is ex- 
actly soluble, and not only ground-state but thermody- 
namics' and dynamicalEP properties have been known. 
S(q, lo) consists only of the two-spinon processes in con- 
trast to the Heisenberg model, and is given by a formula 
analogous to eq. (3.2), as: 

. 1 0(w — uo^ + )Q(lo — tl>£_)0(w u — lo) 



where the boundaries of the continuum are given by 
u>t-(q) = Jq{TT - q)/2, U>e+(q) = J{q - tt)(2tt - q)/2 
and LO u (q) — Jq(2ir — q)/4. 

With these in mind, let us look at the case of a a c . 
Even if a is switched on, S(q,Lo) basically keeps the fea- 
ture of the Heisenberg model within the gapless regime. 
As for the HOS processes, however, their total contribu- 
tion considerably decreases from the value of the Heisen- 
berg model as shown in Fig. 3(a), which shows the sum 
of residues by the HOS processes for small values of q/ir. 
In Fig. 3(b), the maximum values of the above sum be- 
tween < q/ir < 0.5 are plotted as a function of a for 
N = 26. For a ~ 0.2 the HOS contribution becomes 
minimal. 

In view of the conformal field theory, namely as 
for low-energy processes, these two conformal invariant 
models are related to each other as follows. If we cut 
off the exchange coupling longer than the second neigh- 
bor in the HS model, the residual part becomes the J- J' 
model with a = 0.25. The marginally irrelevant opera- 
tors in the discarded long range part cancel out with the 
corresponding ones in the second neighbor coupling of 
0.25 — a c . Thus, it is natural to expect that the behavior 
of S(q, lo) for small lo at a ~ a c resembles that of the 
HS model. Strictly speaking, however, the HOS contri- 
bution never vanishes and the minimum contribution is 
not at a c but ~ 0.18-0.2. This means that the behavior 
of higher-energy processes deviates to some extent from 
the expectation in the infrared limit. 

Such tendency can be seen in Fig. 4, which shows 
S(q,co) for a = 0.2. The shape and the size of the main 
continuum almost coincide with those of the HS model; 
at the upper edge of the continuum, S(q, lo) is more likely 
to have a jump than that in Fig. 1. For more quantita- 
tive discussions on S(q, lo) the introduction of the n-th 
frequency moment defined as, 

/>oo 

K {n) (q) = / dco oo n S{q,co). (3.5) 
Jo 

is convenient. First, we consider static susceptibility 
x{q) = K^ x \q), which is depicted for some values of 
a in Figs. 5(a) and 5(b). 

For the Heisenberg model, the Jiniform susceptibility 
x(0) is given by the Bethe AnsatzJi^P we choose the nor- 
malization as x(0) J = l/(27r 2 ). For a general value oLq, 
diagonalization results have been discussed recently;© 
x(q) is expected to diverges as l/(ir—q) with some log- 
arithmic correction. What we should add here is that 
as in the inset of Fig. 5(a) slight dependence on N is 
observed for small q, and that the extrapolated value for 
q = seems to deviate subtly upward from Griffiths' 
value. This may possibly means the existence of some 
logarithmic singularity in the small-g limit owing to a 
non-fixed-point model. 

As a increases, xio) increases for medium values of 
q. The tendency becomes opposite for the higher-order 
frequency moments, as will be mentioned shortly. As 
discussed above, the curve for a — 0.2, rather than a c , 
coincides well with that of the HS model except for q ~ 0, 
where the long-range part of the coupling probably plays 
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(3.6) 



obtained by integrating eq. (3.4) as, 

x(q)J = \ , 1 . 

^ y/\we-u>e+\ 

In contrast with the Heisenberg case, the leading power 
is linear in q for q — > and x(a)_oc VC 71 " — without 
logarithmic correction for g — > 7T.EZP 

We discuss the cases of a > a c together here. Since 
the system has a finite gap, x(l) vanishes for q — > 0. As 
a sign of it, the abrupt decrease of x(o) can De seen f° r 
large values of a. The value of q at which x{o) drops cor- 
responds to 7r/£, where £ is the spin correlation length. 
For the MG case (a = 0.5), the spin correlation 5j van- 
ishes for i > 2; the anomaly exists at q ~ 7T. By extrap- 
olation we obtain x(7r)J ~ 0.843. 

Leaving the discussion about lO )(g) = 5(g) and 
X^^(g) for §3.3 and Appendix A respectively, here we 
consider K^(q) and K^ 3 \q). As the order n of the fre- 
quency moment becomes higher, the effect of high-energy 
processes is amplified. In Figs. 6(a) and 6(b), the second 
and third frequency moments are plotted, respectively. 
As n increases, the values of the J-J' model for q ~ tt 
are enhanced in comparison with that of the HS model 
given by the integration of eq. (3.5): 



J 2 

1 

~32 
K^(q) 



-^(2tt 
128 v 



q)(7q 2 - 147rg + 67r 2 ) 



(q - 7r) 2 (2g 2 - 47rg + 3tt 2 ) In ( 1 



(3.7) 



J 3 



= ^(2tt- 
768 V 
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(3.8) 



A main cause is probably the HOS contribution of the J- 
J' model outside the TSC, as pointed out for the Heisen- 
berg case.Er The enhancement for small q and large a 
is owing to the considerable intensity beyond the upper 
edge of the main continuum in the gap regime as in Fig. 
3(b). This is one of the topics in §3.3. 

Finally, we point out one more important fact that the 
height of the TSC becomes shorter monotonically with 
increasing a. The highest point of the TSC continuum 
(at q = tt) is wJ for the Heisenberg model, and 7r 2 J/4 
for the HS model, which is comparable with the case of 
a ~ a c . The influence of it appears also in Figs. 6(a) 
and 6(b); as a increases, K^ n \q) for large n generally 
decreases. We return to this point later when we estimate 
the value of J by experimental spectra. 

Having described mainly the resemblance of S(q, u>) for 
a ~ a c to that of the HS model, we close this subsection 
by summarizing the properties of S(q,Lo) for a ^ a c (for 
example, Fig. 4) for later discussions. Basically, the be- 
havior of S(q, lo) remains having common features with 
the Heisenberg model; we emphasize repeatedly [a] the 
intensity of S(q, lo) becomes considerably strong as q ap- 
proaches tt due to the divergence of S(q). [b] The inten- 
sity of S(q, lo) is strongest at the lower edge and a mono- 
tonically decreasing function of lo for every q. As for 



the point [c], although non-two-spinon contribution ex- 
ists, their intensity becomes extremely weak for a ~ a c , 
where the features of S(q,u>) resemble that of the HS 
model. 

3. 3 Excitation spectrum for a > a c 

Since the characteristics of the gap regime are typi- 
cally realized in the MG case (a — 0.5), we take it up 
first in this subsection. In this case we use ty± as the 
starting vectors in the recursion method. In Fig. 7(a) we 
show the residues of C(z) [eq. (2.6)] for N = 24; we plot 
both contributions from ^ + (Q = 0) and (Q = tt) 
together, therefore the spectrum is completely symmet- 
ric with respect to q = 7r/2.Er In Fig. 7(b) S(q,u>) is 
shown for N — 26. We itemize their characteristics in 
comparison with the gapless regime below. 

[1] There is appreciable intensity outside the main con- 
tinuum not only for the small-q region [see also Fig. 3(b)] 
but in the vicinity of q = n. On the contrary, there is 
little high-energy contribution around q = tt/2. 

[2] S(q, lo) is not a monotonically decreasing function 
of lo; the upper edge of the continuum for q > n/2 has 
stronger intensity than the inner area of the continuum. 

[3] The region near the lower edge of the continuum 
has strong weight, like the gapless cases. However, the 
weight of S(q,U)) between q = tt/2 and tt is comparable, 
because the residues become huge around q = tt/2, as 
well as S(q) does not diverge at q = tt. 

[4] An isolated branch appears for 0.45 ^ q/n 0.55. 
As seen in Fig. 7(a) this isolated branch becomes a pure 
single mode for q = n/2 with lo = J independent of 
N. The weak intensity in high energy and the strong 
intensity at the lower edge around q = tt/2 mentioned 
in [1] and [3] respectively are attributed to this isolated 
branch. 

These features are in sharp contrast with the ones [a]- 
[c] for a < a c . 

Let us look at the spectra more closely. The structure 
of the continuum is more complicated than the TSC in 
the gapless regime; near the upper edge it is not easy to 
trace each branch, while a series of dominant branches 
can be distinguished in the lower half of the continuum. 
In particular, near the lower edge, a number of branches 
with comparable weight seem to gather. The lowest 
branch does not have the strongest intensity. 

To check the dependence on N, we show the residues 
of the four lowest dominant branches in Figs. 8(a) and 
8(b). For q/ir 0.55 the weight of every branch depends 
severely on system size. Although the branches 1' and 2' 
monotonically decrease and the branch 4' monotonically 
increases with increasing N, the branch 3' shows compli- 
cated behavior. Positions of the corresponding poles are 
shown in Fig. 8(c). The higher the energy of the branch 
is, the severer the dependence is. As shown in the inset 
of Fig. 8(c), which shows the pole positions versus 1/A 2 
at q = 7T, every branch seems to converge to the lower 
edge. 

It is clear that the spectrum forms a continuum for 
q/ir J; 0.55, because both the residues and the poles are 
considerably dependent on system size. On the other 
hand, as in the Figs. 8(a) and 8(b), for q/ir ^ 0.55 the 
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weight of the branches 1' and 2' is enhanced rapidly as q 
approaches tt/2 without dependence on N. Inversely, the 
branches 3' and 4' seem to vanish. As for the pole posi- 
tions of the branches 1' and 2', there is little dependence 
on N [Fig. 8(c)]. Consequently, according to the discus- 
sion in §3.1, an isolated branch with weight more than 
90% seems to occur in this region. In fact, at q — tt/2 it 
is proven that the spectrum becomes a pure single mode 
as shown in Appendix B. Incidentally, a similar spectrum 
was obtained by using a variation method;© in addition 
to a continuum, a bound state below the continuum was 
derived for 0.36 £ q/ir £ 0.64. 

Next, we touch on the point [2]. In Fig. 9 the residues 
at q = tt are plotted versus ui. Near the lower edge, 
there exists a large peak like the gapless cases. However, 
it is not clear in this case whether this peak diverges, and 
whether it is situated exactly at the lower edge. Anyway, 
the weight in the continuum once has a broad minimum 
(oj/J ~ 1.5-1.8), then another peak appears near the 
upper edge (uj/J ~ 2). This peak seems less dependent 
on N, and has tails in both sides of ui. Similar tendency 
can be seen for other values of q. This profile is in sharp 
contrast with those in the gapless region. 

Having focused on the MG case, here we mention the 
cases of smaller values of a. In these cases, the ground 
state is not twofold degenerate except for the thermo- 
dynamic limit, thus we treat excitations only from the 
nondegenerate ground state for finite N . In the thermo- 
dynamic limit, the two ground states with Q = and tt 
will give contributions mutually symmetric with respect 
to q = 7r/2 just like the MG case. 

Figures 10(a)-10(c) show S(q,uj) for TV = 26 and a = 
0.35, 0.4 and 0.45, respectively. In these figures we can 
recognize that the characteristics of the MG model [l]-[4] 
still remain. Comparing these with Fig. 4 (a = 0.2) and 
Fig. 7(b) (a = 0.5), we find that the spectrum gradually 
changes from the gapless feature to the one of the MG 
case. 

Finally, for later discussions, we represent some static 
quantities; Fig. 11(a) shows the ground-state energy, and 
Fig. 11(b) the lowest excitation at q = tt/2 and the gap 
A. These values are extrapolated from the finite-size 
data by the polynomial fit of second and fourth orders. 
There is no critical behavior at a = a c . 

The static, structure factor S(glhoth for the Heisen- 
berg modeled and the HS modeKtS 1 is logarithmically 
divergent at q = tt. Accordingly, the situation is basi- 
cally the same for the gapless regime of the J- J' model. 
On the other hand, as a increases beyond a c , the di- 
vergence at q — tt is suppressed. For the MG point 
S(q) = (1 — cosg)/4, namely S(q) is not only conver- 
gent but upward convex near q = tt. In Fig. 12 S(q) for 
large a is shown. Even for a — 0.4, S(q) is considerably 
enhanced near q = tt. Furthermore, it increases with in- 
creasing system size. On the other hand, for a — 0.45 
the dependence on N of S(tt) is in the direction of sup- 
pression. 

Anyway, concerning the excitation in the regime of 
a > a c , we are still far from understanding in the stan- 
dard of the Heisenberg model. Further studies are nec- 
essary. 



§4. Comparison with Spin-Peierls Compounds 

Based on the discussion in the preceding section, we 
consider various aspects of the SP transition in CuGeOa 
and a'-NaV 2 05. 

4-1 Complete spectrum of inelastic neutron scattering 

Recently, a complete spin excitation spectrum of 
CuGeOs was obtained by the inelastic neutron scattering 
by Arai et alS3 The gap and the size of the continuum 
in this spectrum |ajje_cansistent with the previous neu- 
tron experiments.E2Hl3'E2P At a glance on the spectrum 
for T = 10K (Fig. 1 of ref. 19), it looks like the TSC of 
the Heisenberg model. However, as the authors pointed 
out, there arc points which are incompatible with the 
Heisenberg model. We itemize such points below includ- 
ing the authors' words. 

[0'] The shape of the lower and the upper edges of the 
continuum is wedge- like and does not have the roundness 
characteristic of sine (a = 0) and quadratic (a ~ a c ) 
curves. Concerning this point, compare Figs. 18(a) and 
19(a), for instance. 

[1'] There is appreciable intensity above the contin- 
uum, especially for q ~ tt. 

[2'] There is "rampart" or ridge of scattering sur- 
rounding a valley in the spin continuum. 

[3'] There exists a peak in intensity on the lower 
boundary of the continuum at zone boundaries (0,0,7r/2) 
and (0,0,3tt/2). 

The points [l']-[3'] completely correspond to what we 
have summarized as [l]-[3] for a > a c in §3.3, respec- 
tively. In this viewpoint, let us compare Fig. 1 of ref. 19 
with Figs. 7(b) and 10(a)-10(c) in this paper, before tak- 
ing account of the effect of coupling alternation. Among 
the above four panels Fig. 10(c) (a = 0.45) seems to bear 
the best resemblance. 

As for the point [1'], the intensity of every pole over the 
upper edge of the main continuum (q ~ tt) in Fig. 10(c) 
corresponds to the one in Fig. 1 of ref. 19. For example, 
a horn-like short branch jutting out upward from q = tt 
just above the main continuum (uj/J ~ 2.1) is distinctly 
reproduced. Incidentally, because the spectrum should 
be symmetric with respect to q = tt, the asymmetric 
intensity in the experimental spectrum, especially for q > 
n, is due to the contamination by phonons etc., as the 
authors mentioned. 

Concerning the point [2'], the strong intensity at the 
upper edge and the weak intensity in the upper half of 
the continuum for q ~ 27r/3 - tt are quite similar between 
the two figures. Recall the profile in Fig. 9 for the MG 
case. 

As for the point [3'], the experimental intensity near 
the lower edge is almost the same between q = tt/2 and 
tt. Furthermore, the maximal points seem to be situated 
slightly above the exact lower edge, which we have sug- 
gested in §3.3, although it is not clear whether this is due 
to its intrinsic property or thermal broadening. 

For smaller values of a [Figs. 10(a) and 10(b)], the 
tendencies deviate from the experimental spectrum. For 
example, in Fig. 10(a) for a = 0.35, the strongest in- 
tensity, as to the lower edge, is clearly at q = tt. If 
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a ~ a c , the intensity over the continuum considerably 
attenuates, as discussed in §3.2. Consequently, the value 
of a should be fairly larger than a c . 

More quantitative fit is possible; in Fig. 13 we com- 
pare integrated amplitude S(q) (Fig. 12) with the neu- 
tron data, which are taken from Fig. 2 of ref. 19. In 
the experimental results, there is neither divergence nor 
abrupt increase for q ~ tt characteristic of smaller val- 
ues of a. It reads a ~ 0.45 for T=10K and a ~ 0.5 
for T = 20K (T SP = 14K), which is consistent with the 
above direct comparison of S(q,ui). This indicates that 
the value of a hardly changes through Tsp, or reduces 
slightly for T < Tgp if any. Incidentally, there is a ten- 
dency of the experimental data to be rather smaller than 
the diagonalization result for q ~ n. This is probably due 
to two facts. One is that in obtaining S(q) by integrating 
experimental S(q,u>), the higher- frequency contribution 
than 32meV (the top of the continuum) is cut off to 
exclude the contribution by phonon. Actually, there re- 
mains magnetic contribution of several percent, as men- 
tioned above. The other is that the introduction of 6 
somewhat reduces S(q) at q = tt. In this connection, if a 
smaller value of a is assumed, a considerable large value 
of 6 is needed to reproduce the experiment, for exam- 
ple, 5 ^ 0.03 even for a = 0.4; such a parameter set is 
incompatible with the value of A, as we will see later. 

Thus, we have found that most features of this neutron 
spectrum are explained by the J-J' model with a = 0.4- 
0.45, even without introducing S. 

Next, let us check whether or not the above aspects 
are also realized by introducing 6 without a or with 
smaller values of a. Figures 14(a)=14(c) show S(q,ui)'s 
for 5 = 0.05 and three values of ao To emphasize the 
effect of 8, a fairly large value of S is used for CuGeOs 
(8 is of order 0.001). In these figures some new features 
can be seen. For example, as will be discussed in §4.3, 
there appears an opening of intensity between the lowest 
branch and the above continuum, especially for large a. 
However, for the case of a < a c [Figs. 14(a) and 14(b)] 
the characteristics [l]-[3] for a > a c discussed in §3.3 do 
not appear even if 5 is added, while these characteris- 
tics are preserved for a > a c if 5 is added, as in Fig. 
14(c). In other woxds, the characteristics [l]-[3] appear 
only when a > a c .c3) Furthermore, these characteristics 
do not become manifest unless a is fairly larger than a c . 
For smaller values of 5, the characteristics [l]-[3] do not 
appear for a < a c , of course. Thus, these facts exclude 
at least the possibility of a < a-, which is contrary to 
the discussion of Castilla et alt3 We will point out a 
mistake in their assumption in §4.3. 

4-2 Estimation of J and 6 

First, we discuss the value of J. The first estimation 
of J was given by a neutron experiment as J ~ 120K,EP 
in which the maximal point of the lowest branch in the 
spectrum was compared with Jtt/2, which is obtained 
from the des Cloizeaux-Pearson curve [the first equation 
of eq. (3.1)] for the Heisenberg model. However, the 
experimental spectrum has to be compared with that of 
eq. (2.1). As we have already shown in Fig. 11(b), the 
size of the continuum decreases with increasing a. The 



maximal value of the lower edge (q — 7r/2) decreases from 
Jtt/2 to J, as a changes from to 0.5. Since a = 0.4- 
0.45, J is modified into a larger value by ^50%; thus we 
have J - 180K. 

In the above estimation of J we have neglected the 
effect of 6; let us check it below. In Fig. 15 the maximal 
values of the lower edge are plotted as a function of S 
for some values of a. In contrast with the effect of a, 
the spectrum size weakly depends on S. Although it 
abruptly increases for small S and large a, the increment 
is at largest a few percent, because S is extremely small. 
Consequently, the coupling alternation scarcely affects 
the estimation of J. 

Now, we refer to the estimations of J by other means. 



Many authors tried to evs 



uate a and J by comparing 



with x(T) by experiment,EP which shows a broad peak 
around 56K. Among them, Riera and Dobry obtained 
a = 0.3£Land J = 160K, by fitting the maximal point 
of %(T)0 Since this fit was quite good for T > 56K, 
following studies often adopted these values. 

Using Faraday rotation, Nojiri et alx3 estimated J 
at 183K; they compared the field at which magnetiza- 
tion was saturated with the formula for the Heisenberg 
model: g^H s — 2J. Although in this case the corre- 
sponding value for eq. (2.1) should be used for accuracy, 
it is convenient_to use, instead, the known formula for 
the J-J' modem 



gfj, B H s /J = 



2 

l + 2a+l/8a 



< a < 0.25 
a > 0.25 



(4.1) 



Unless the variation by 5 is singular at S = 0, this esti- 
mation is not poor because 8 is very small. As a result, 
the value is slightly modified as J = 178K, 173K, 168K 
for a = 0.35, 0.4, 0.45, respectively. 

Thus, the estimation of J by the neutron experiments, 
J ~ 180K, is consistent with other estimations. By us- 
ing a — 0.45 (0.4) and J = 180K, the gap energy with- 
out coupling alternation is evaluated at 19K (5K), which 
must be a value too small to make a gap by itself above 

Next, we estimate the coupling alternation parameter 
6. The excitation gap was estimated by some experi- 
ments. In particular, the neutron scattering experiment 
can directly observe the gap as a function of temperature 
A(T), and concluded A(0) = 2.1meV (24K)Jj3> which is 
consistent with other experiments like. NMR.pa' specific 
heatjO) magnetic susceptibility x(T)jS etc. 

In Figs. 16(a) and 16(b), the gap A versus 6 obtained 
by the exact diagonalization for eq. (2.1) is depicted. In 
Fig. 16(b), the dash-dotted (dashed) line corresponds to 
A = 24K and J = 180K (160K). If we assume a = 0.45 
(0.4) and J = 180K, 6 becomes as tiny as 0.001 (0.005). 
This value is muph | smaller than the previous estimations: 
S = 0.012-0.014HE3) This difference is originated in the 
fact that A is an abruptly changing function of a for 
~ 0.5 and of 5 for ~ 0. 

Experimentally, the atomic displacement due to the 
lattice dimerizatiqn was measured for CuGeOa by neu- 
tron diffraction.EP The displacement ratio Si is very 
small; 8i ~ 0.001 both for Cu along the c axis parallel 
to the ID chain and for 0(2) along the axes perpendic- 
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ular to the chain. i-This value is too small to detect in 
some other means. O However, it is still unknown how 
this small displacement is connected to a small exchange 
coupling. 

In CuGeC>3 , CuOg octahedra link mutually with edge 
sharing along the c axis; the angle between Cu-O-Cu 
is 98°. The exchange coupling, which is determined by 
the overlap of the involved localized orbitals, tends to 
be weakly ferromagnetic in the case oLthe right angle, 
as Kanamori classified many years ago-E^P On the other 
hand, recently Geertsma and Khomskii argued that for 
the realization of antiferromagnetic coupling the side- 
group effect of Ge plays an important roleJiZP Thus, we 
cannot assert the relation 8g_ oc 8 necessarily holds. 

By the way, from the above discussion one can ex- 
pect relatively large long-range exchange couplings sur- 
vive in edge-sharing materials. Hence, we have checked 
the third-neighbor's effect for various strength of 
(> 0) and for some values of a and S. The common 
aspect among the resultant spectra is that the top of 
the lower edge of the continuum becomes broad and flat 
around q — tt/2 with the spectral size almost unchanged. 
This aspect is quite different from the characteristics of 
CuGeOs. Thus, the third-neighbor exchange is found 
negligible. 

Now, we refer to the second inorganic SP compound 
a'-NaV205j!3' in which VO5 pyramids link mutually 
with corner sharing in the ID direction. Since x(T) for 
T > Tgp is well fitted by the Bonner-Fisher curveH with 
J =560K, the frustration by J' seems small in this com- 
pound. The gap is estimated at A = 98K by NMR.EP 
Using these values or another estimation of J = 441K 
and A = 85K,E$ we have^I = 0.043-0.048 from the 
curve of a = in Fig. 16(a) .Et Since the relation 81 oc 8 
holds in this case, the displacement due to dimerization 
is probably much larger. In fact, a split in the Na- 
NMR spin-echo intensity corresponding to the chain di- 



qcn observed below Tsp,E_F in contrast with 



rection ha 
CuGe0 3 .E 

Before closing this subsection, let us look at the SP 
transition in the light of energy. We depict the total 
magnetic energy E/J for 8 = and 8 ^ in Figs. 11(a) 
and 17, respectively. As temperature is lowered, the lat- 
tice will be distorted so as to minimize E/J. Thus, the 
displacement at T = Tgp will tend to enhance S as well 
as J and lower a. This tendency agrees with the one 
seen in Fig. 13, although the thermal fluctuation may be 
also involved in this case. 



4-3 Further discussions 

In this subsection, we will consider a couple of issues 
in relation to the experiments. ._. 

First, we take up A(T). Harris et allEy discussed this 
issue by combining three materials: (1) The dependence 
on temperature of the intensity of a superlattice peak: 
1 oc [Tsp — T] 2/3 . The work by neutron diffraction 
f and X-ray diffractiontS concluded (3 ~ 0.33 in 
agreement. (2) The spontaneous contraction along the b 
axis with the Ginzburg-Landau theory: I oc 8f. (3) The 
Cross-Fisher theory: A oc <5 2 / 3 . Thus, they obtained 
A oc [Tsp — T] 2a " with a — 1/3, tacitly assuming S oc Sg. 



Later, Castilla et aL0 excluded the possibility of 
a > a c , based on the above discussion and the 
renormalization-group theory. They assumed A oc <5 2 / 3 
with negligible logarithmic correction throughout the 
gapless regime, while A oc 8 for a > a c . However, this 
assumption is false. As in Fig. 16(a), A is nearly pro- 
portional to S for small 8 and a — 0. The leading power 
<5 2 / 3 is smeared out due to the severe logarithmic correc- 
tion, even for large values of 8. This aspect had. 
already pointed out by the diagonalization studies .t 
On the contrary, the leading power (5 2 / 3 was clearly ob- 
served for a = a c in the density matrix renormalization 
group calculations by Chitra et al. ,Ef because the model 
becomes conformal invariant and there is no logarithmic 
correction, as mentioned in §3.2. Our data for a — 0.2 is 
consistently fitted with a power of ~ 0.71. As a becomes 
larger, the seeming power further decreases with a finite 
gap at 8 = and approaches 0.5 near the MG point. 
These results are rather contrary to the assumption by 
Castilla et al. Consequently, there is no theoretical rea- 
son that a has to be smaller than a c . 

Actually, the jalue a in the power of A(T) was esti- 
mated at ~ 0.15 jE^f which is much smaller than the seem- 
ing power in Fig. 16 of ~ 0.26 for a = 0.45, as well as the 
prediction of the Cross-Fisher theory of 1/3. However, 
the discussion of power itself is meaningless -for the gap, 
which has serious logarithmic modification.E3' Further- 
more, it is not certain whether the relation oc 8 holds 
for edge-sharing materials like CuGeOa , as discussed be- 
fore. On the other hand, A(T) for o/-Na.V205 also shows 
a similar abrupt increase at T = TgpEy These may in- 
dicate that A(T) is independent of the details of the 
exchange coupling. 

Second, we consider the double gap xeported in a re- 
cent neutron experiment for CuGeOsH According to 
it, gap- like weak intensity was observed above the low- 
est peak at 2.1meV for q = ir; this "second gap" had a 
similar width with the first one. 

Let us return to Figs. 14(a)-14(c) for 8 — 0.05. Com- 
paring these figures, we recognize that the aspects are 
different between Figs. 14(a) [or 14(b)] and 14(c). In the 
former case (a < a c ) the lowest branch for 8 = [Figs. 1 
and 4] splits into two and the opening between them be- 
comes wider with increasing a; at q = n, however, these 
two branches merge and the second gap seems to vanish. 
Meanwhile, in the latter case (a > a c ) the second gap is 
open for all the range of q. Thus, this issue is highly de- 
pendent, at least quantitatively, on the value of a and 8. 
Here, we focus on the plausible parameters for CuGeOa 
and a'-NaV205, leaving more general discussions for fu- 
ture publication. 

To begin with, we consider the case of a — 0.45, 
namely CuGe03. Shown in Figs. 18(a) and 18(b) is the 
dependence on N of the pole position and the residue re- 
spectively corresponding to the two lowest branches for 
S = 0.01. We use a somewhat larger, but sufficiently 
small, value of 8 than our estimation so as to find its 
effect. 

First, we discuss the branch 1'. Both the pole posi- 
tion and the residue of the branch 1' scarcely depend on 
system size except for q ~ or n. For q = n, relatively 
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reliable extrapolation for both quantities is possible, as 
shown in the inset of Fig. 18(a) and by an arrow in Fig. 
18(b). Moreover, since the pole position of the branch 2' 
is always away from the branch 1', as will be discussed 
below, the branch 1' forms an isolated mode. 

Next, as for the branch 2' both the pole position and 
the residue are almost independent of system size for 
medium values of q, although the behavior of the residue 
is complex. For q = tt, the pole position converges to 
a value which is higher than that of the branch 1', as 
shown in the inset of Fig. 18(a). On the other hand, we 
are not certain that the residue of the branch 2' at q = tt 
has finite weight for N — > oo. Concerning the branch 3' 
(not shown in the figures) , both the pole position and the 
residue depend on N, but we cannot elicit definite an- 
swers from our data whether its pole position converges 
to that of the branch 2'. Consequently, the branch 2' 
probably forms another isolated branch, although there 
remains the possibility of the lower edge of the contin- 
uum. The branch 3' belongs to the continuum. 

Anyway, the opening between the branches 1' and 2' 
exists for every value of q. This opening is thin for 
medium values of q and relatively wide at q — tt. Nev- 
ertheless, the width of this second gap at q = tt is fairly 
smaller than that of the first gap, as shown in the inset 
of Fig. 18(a). Incidentally, if the branch 2' is an isolated 
branch at q = tt too, the third gap will also exists. 

We have confirmed this situation remains the same 
for some similar values of the parameters. On the other 
hand, for larger values of 5 (for example a — 0.45 and 5 = 
0.1), the branches LLand 2' form clear isolated modes, 
as studied formerly.EZP Even in this case, however, the 
width of the second gap is found still smaller than the 
first gap. This result is in sharp contrast with an RPA 
calculation)!^ in which these two gaps have the same 
width. 



Now, we turn to t 
a'-NaV2 05 in mind. 



case of a — and S = 0.05 with 
3* In Fig. 19(a) we show the pole 
positions of the lowest three branches with dominant in- 
tensity. The branches 1' and 2' are almost independent 
of system size for all the range of q. A conspicuous point 
is that these two branches merge into a single pole at 
q = tt. On the other hand, the poles belonging to the 
branch 3' have severe dependence on system size. More- 
over, we find by the extrapolation of the poles at q = tt 
that the branch 3' converges to the branch 2' [the inset 
of Fig. 19(a)]. Shown in Fig. 19(b) are the corresponding 
residues. As for the branch 1', the dependence on N is 
seen for q/ir ^ 0.6; the branch 2' also has large depen- 
dence on N. Consequently, we can interpret the above 
as follows: the branch 1' is isolated except for q ~ tt, and 
the branch 2' is the lower edge of the continuum. Thus 
the second gap vanishes for q = tt, while an isolated 
branch appears appreciably below the continuum for q 
away from q = tt. Furthermore in the opening between 
them there is weak intensity of the higher-order processes 
for q/ir ~ 0.5 ± 0.2. The situation is quite different for 
larger values of 8 ( ^5 0.15), for which there appears an 
obvious second gap even at q — tt. However, the second 
branch belongs to a continuum unlike the cases of large 
a. 



Summarizing, the appearance of a second gap or iso- 
lated branches is dependent on the values of a and 5. The 
second gap, which is narrower than the first gap, may be 
found at q = tt for CuGeOa. On the contrary, an isolated 
branch will be observed for q ^ tt for a'-NaV205. 

Finally, we mention a connection with the neutron ex- 
periments under high pressure.Ea-' T$p increases linearly 
with increasing pressure The lattice constant along 
the chain (c axis) becomes slightly longer, while the dis- 
placement of oxygen in the ab plane is rather large. The 
gap energy becomes twice larger (4.1meV) under 1.8GPa 
at 5K. On the other hand, the dimerization for Cu atoms 
decreases under high pressure. 

One possible interpretation of the above experiments, 
which seem mutually contradictory at a glance, within 
the model of eq. (2.1) is that a is enhanced almost up 
to the MG point by applying pressure. If so, the gap 
becomes about twice larger even with reduced 6 as in 
Fig. 16(b). Since the abrupt enhancement of A is lim- 
ited to the vicinity of the MG point as in Fig. 11(b), this 
fact may be an evidence of the strong frustration. Any- 
way, more detailed experiments and reliable theoretical 
knowledge on the relationship between the lattice struc- 
ture and the model parameters are needed to settle this 



§5. Summary 

Using the exact diagonalization and the recursion 
methods, we have investigated dynamical as well as static 
properties of the ID J-J' model with coupling alterna- 
tion, and compared with experiments of inorganic spin- 
Peierls compounds, CuGeOa and a'-NaX^Os. First, 
we have discussed the dynamical properties of the J- 
J' model without coupling alternation. The main points 
are: 

(1) For a < a c , the characteristics ([a]-[c] in §3.1 and 
§3.2) of the excitation spectrum are basically the same 
with those of the Heisenberg model. 

(2) For a ~ a c , the spectrum resembles that of the 
Haldane-Shastry model. 

(3) As a further increases beyond a c , there ap- 
pear qualitatively different features from the Heisenberg 
model, summarized as [l]-[4] in §3.3. 

Next, based on these results, we have studied a variety 
of features of the above SP compounds, particularly in 
view of the inelastic neutron scattering experiments. We 
itemize noticeable points below: 

(4) Most characteristics of the complete neutron spec- 
trum at 10K are quantitatively reproduced by the ID 
J-J' model with a = 0.4-0.45. Some experimental as- 
pects cannot be explained, unless a is at least somewhat 
larger than a c . 

(5) Using a — 0.45 (0.4), the strength of exchange 
coupling is evaluated at J ~ 180K, consistent with the 
estimations by other means. The dimerization parameter 
of the exchange coupling 5 is found extremely small as 
- 0.001 (0.005). 

(6) We have investigated the isolated branch appear- 
ing below the continuum or the double-gap structure for 
CuGeOa and a'-NaV2 05. The features of the second gap 
will be different between these two compounds. 
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(7) Pressure effect on CuGeC>3 is possibly interpreted 
as the increase of the frustration. 

Our estimation oLa is somewhat larger than the ones 
by x(T), oe ~ 0.35 J2J^ Some causes of this discrepancy 
can be thought of. For example, evaluations are made 
at different temperatures between the two means. Our 
estimation does not allow for the temperature effect. 
There is quantitative ambiguity in neutron experiments 
at the present stage. Leaving aside such a quantitative 
point, we would like to emphasize that the effect of the 
large frustration beyond a c is essential for various unique 
properties of CuGe03. 

There remain many significant issues left for further 
studies. For example, although we treat a pure ID model 
because the SP transition itself proves good one dimen- 
sionality, a relatively large magnetic dispersion perpen- 
dicular to the ID chain was observed in CuGe03.tHf How 
do we interpret it? What effect does the large frustration 
have on the elemental substitutions? Furthermore, pll£j- 
non softening has never been discovered in CuGe03,EP 
which has always been a driving force of the organic SP 
transition so far. Does the large frustration-shoulder such 
a role, as conjectured by Buchner et al. IcB To this end, 
the accumulation of experimental data on a'-NaV2 05 
will be important. 
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Appendix A: Single-Mode Approximation 

In this Appendix, we summarize the— properties of 
the single-mode approximation (SMA^'El in the J- J' 
model. As a kind of variational trial state for a spin- 
excited state with momentum q we consider 

m = s z q \o), (at) 

where |0) is the ground state — in this case exact. The 
excitation energy with respect to tp q is written as 

<0| [S*_ q ,[H,S*]}/2\0) 



E„ — E 



E 



(^ q ) ^ (Q\si q s»\o) 

KW(q) _ JduojS{q,oj) 

~ fdu;S(q,u;) = 



where K^(q) is the first frequency moment. As seen in 
the last expression of eq.(A.2), iv q can be regarded as a 
kind of average of w [the distribution function is S(q, u>)], 
besides as a variational excitation energy. Hence the 
SMA becomes good when the distribution of S(q, u>) con- 
centrates on a small range of u; in particular if the weight 
of S(q, to) is restricted to a certain sole value, the SMA 
becomes exact. We will prove that a pure single mode 
is realized for the exactly soluble cases at q — ir/2 in 
the next Appendix. Inversely, when the distribution of 
S(q) spreads over a wide range of u> evenly, the SMA is 
poor as the approximation of the lowest mode. In the 
meantime, since ui q is an average, it is an upper bound 
for the lowest spin excited mode. Therefore, one can use 
the SMA for a proof of the gaplessness. 

For the Hamiltonian eq. (1.1), the first frequency mo- 
ment is calculated as 

K^(q) = -2J[(l-cosg)5i + a(l -cos2g)Sy , (A-3) 

where Si is the two point correlation function with re- 
spect to the ground state: 

1 



Si ~ N ^2^jSj +i ) 



(A-4) 



Thus, by using the averages of the exact diagonalization, 
the SMA can be estimated. In Fig. 20 we depict the 
results for the J- J' model, thus obtained. 

For a — and small q, the SMA curve is slightly higher 
than both of the two curves surrounding the TSC: this 
is due to the HOS processes as mentioned in §3.1 .EZP For 
a = 0.2 the result is close to that of the HS model, for 
which the SMA is estimated througliJthe known analyt- 
ical formulae of JfW(g)EP and S(q) B> 



J 



16 



(27r-«) + ±( g -7r) a ln(l-|), (A-5) 



S(q) = - 



(A-6) 



This tendency has been discussed in §3.2. 

For the MG case, (q) = J(l - cos q) /4 and S(q) = 
(1 — cosg)/4, therefore 0j q = J independent of the value 
of q. As will be shown in the next Appendix, the spec- 
trum becomes purely single at q = ir/2; the SMA be- 
comes exact then. Similarly, the SMA is relatively good 
for a — 0.45 and q ~ 7r/2, though not exact. On the 
other hand, for q = or tt the SMA is poor as the low- 
est excitation. In the viewpoint of an average, however, 
the SMA indicates that there exists sizable contribution 
above u> = J to balance the lower excitation modes, as 
seen in Fig. 7(a). 

To discuss the existence of the gap quantitatively in 
the SMA, larger system sizes are desired. 



Appendix B: Proof of a Single Mode for Ex- 
actly Soluble Cases 

In this Appendix, we describe a proof a la Majum- 
(A-2) dar an d GhoshtP that S(q,u>) at q = tt/2 in the exactly 
soluble cases (2a+S = 1) is a pure single mode with a ex- 
citation energy independent of system size. The case of 
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(5 = has been already discussed ^ and a paore compact 
proof was given of ipe being an eigenstateES' 

If Sgl^o) becomes an eigenstate of H, the surviving 
contribution to the sum of eq. (2.4) is given only by this 
eigenstate, due to the orthogonality of the eigenstates. 
Therefore, what we have to do for the proof is to show 
Sg\^o) is an eigenstate of H, and then to estimate its 
eigenvalue. 

The proof that tp c = S^L^I^o) becomes an eigen- 
state of H is obtained by an -analogous treatment with 
the proof of the ground stateilP Henceforth, we consider 
systems of N = 41 (I: integer). Before going to the 
proof, we summarize notations and relations. Like a sin- 
glet pair of eq. (2.3), we denote a triplet pair by 



{i,j} = a(i)/3(j) + P(i)<x(j)- 
Useful relations for this proof are as follows!} 1 

^(l-iS i -S m )[£,m] = 2[£,m], 
l(l-ASf8 m ){t,m} = 0, 



- (1 - iSg ■ S m ) [k, £} [m, n] = [£, m] [n, k] 



(B-l 

(B-2 
(B-3 
(B-4 



(1 
(1 

[k, £] [m, n 



4,S e 



S m ) [k,£]{m,n} 

S m ){k,£}[m,n] = -[£,m]{n,k}, (B-5 
I- [k,n] [£, m] + [k,m] [n, £} = 0, (B-6 

[k, £]{m, n} + [£, m]{n, k} + [m, n]{k, £} + [n, k]{£, m} = 

(B-7 

First, we consider the case of 6 = 0. To make use of 
the above relations, we rewrite the Hamiltonian eq. (1.1 

as, 



with H = Hi 



H=jNJ(l + a)-^JH, 
f H2 and, 

N -, 



3+1 ) 



H 2 



N 1 

3=1 



AS* ■ S 



3+2 J 



(B-8 



(B-9 



(B-10 



As a ground-state eigenfunction, we take a form of eq. 
(2.2); substitution by ip2 makes no difference. Then, the 
function of our concern is given as, 



h = S'-\il>i) 



1 



N 



5Z ex p 



7T 



= -7= {lS{ 



5, 2 



1 



2VN 



S| + --0 [1,2] 
(1 + i)(ri - r 3 + t 5 - t 7 H 



TN-3 ~ Tjv-1 



■ ■■[N-1,N] 
(B-ll) 



Tj = [1,2] • • • \j-2,j-l]{j,j+l}\j+2,j+3] ■ -'[N-1,N], 

(B-12) 

with odd j. 

First of all, let us take up Hti, for example. Applying 
Hi and H2 to t\ and r 3 and using eqs. (B.2)-(B.5), we 
obtain 

Hxr x = (N - 2) n - [2, 3]{4, 1}[5, 6] • • • [N - 1, N] 
+{l,2}[4,5][6,3]---[7V-l,iV] + --- 
+{1, 2}[3, 4] • • • [JV - 2, N - 1][N, N - 3] 
— [N, 1]{N - 1, 2} [3, 4] • • • [N - 3, N - 2], (B-13) 

7Yir 3 = (iV-2)r 3 -[2,3]{4,l}[5,6]---[iV-l,iV] 

-[1,2][4,5]{6,3}[7,8]---[JV-1,JV] 

+ [1, 2]{3, 4}[6, 7] [8, 5] [9, 10] • • • [JV - 1, N] + ■ ■ ■ 

+ [1,2]{3,4}[5,6] ■ • • [JV — 2, JV — 1][JV, JV — 3], (B-14) 

H 2 T 1 /a= - [1, 3]{4, 2}[5, 6] • • • [N - 1, JV] 
+ [2,4]{3,1}[5,6]---[JV-1,JV] 

-2{1, 2} [3, 5] [6, 4] [7, 8] • • • [JV - 1, JV] 

-2{1, 2} [3, 4] • • • [JV - 3, N - 1] [JV, JV — 2] 



+ [N- l,l]{JV,2}[3,4]---[JV-3,JV-2] 
- [JV, 2] { JV - 1, 1} [3, 4] • • • [JV - 3, JV - 2] , 

W 2 T 3 /a = [l,3]{4,2}[5,6]---[JV-l,JV] 

-[2,4]{3,1}[5,6]---[JV-1,JV] 

-[1,2][3,5]{6,4}[7,8]---[JV-1,JV] 

+ [1,2][4,6]{5,3}[7,8]---[JV-1,JV] 

-2[1, 2]{3, 4} [5, 7] [8, 6] [9, 10] • • • [JV — 3, JV — 2] — 

-2[JV - 1, 1][2, JV]{3, 4} [5, 6] • • • [JV - 3, JV - 2]. 



(B-15) 



(B-16) 



To begin with, we consider terms appearing in H2(j\ — 
r 3 ). The first two terms in eqs. (B. 15) and (B.16), 

-2 ([1, 3]{4, 2} + [4, 2]{1, 3}) [5, 6] ■ • • [JV — 1, JV], 

are transformed into 

2 ([3, 4]{2, 1} + [2, 1]{3, 4}) [5, 6] ■ • ■ [JV - 1, JV], (B-17) 

by using eq. (B.7). We assign (or should say return) the 
first term of eq. (B.17) to H2T1 instead of the first two 
terms of eq. (B.15), and the second term of eq. (B.17) to 
H2T3 instead of the first two terms of eq. (B.16). Simi- 
larly, we can substitute 

2[JV-l,JV]{l,2}[3,4]..-[JV-3,JV-2], 

for the last two terms of eq. (B.15) by considering 
H2(—t~n-i + ri). Next, we apply the relation eq. (B.6) 
to the two singlet pairs on which H.2 has been operated 
(namely pairs of irregular order) in residual terms of eq. 
(B.15). This substitution generates diagonal terms (ti) 
and the same kind of terms appearing in H\T\ [from the 
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third term to the last but one in eq. (B.13)]. The residual 
off-diagonal terms, namely the second and the last ones 
in eq. (B.13) cancel out with the corresponding terms in 
eq. (B.14) and HiTn-i, respectively. Consequently, eqs. 
(B.13) and (B.15) are modified as 

(Hi + H 2 )n = (N- 2)n + aNn 

+ (1 - 2a) ({1, 2}[4, 5] [6, 3] [7, 8] • • • [N - 1, N] + ■ ■ ■ 

+{1, 2} [3, 4] • • • [N - 2, N - 1][N, N - 3]) (B-18) 



If we put a = 1/2, the off-diagonal part in eq. (B.18) 
vanishes, resulting in 



2 Ti. 



(B-19) 



As for T3, T5, • • •, the same form with eq. (B.19) can be 
derived in the same manner. Then, by considering eqs. 
(B.ll) and (B.8), the following eigenvalue equation can 
be obtained for the Majumdar- Ghosh model, 

HVe= (-|jV+l)j^ e . (B-20) 

Since the ground-state energy is Eq = —37V J/8, the 
eigenvalue of tp e = S^ 2 ipi is above Eq by J, indepen- 
dent of the system size N. 

This proof can be extended to the case of the exactly- 
soluble dimerized J- J' model eq. (2.1) with 2a + 8 = 1. 
Note that ip2 cannot be chosen in eq. (B.ll) here. Then, 
the eigenvalue equation corresponding to eq. (B.20) be- 



-~(l-a)N + 2(l-a) 



Jlpe- 



(B-21) 



Since, the first term is the ground-state energy in this 
case,fj the energy increment is 2(1 — a) J. 
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Fig. 1. S(q, u>) of the Heisenberg model for TV = 26. The poles are 
situated at the center of the spheres. The intensity of each pole is 
proportional to the volume of the sphere to emphasize the poles 
with weak intensity. The boundaries of the two-spinon contin- 
uum, u>e(q) and uj u (q), are shown by solid lines. The lower edge 
of the higher-order-spinon continuum is indicated by a dashed 
line (TV = 26), which is strongly dependent on TV. The lowest 
quintet level for each q (TV = 14) is marked with solid diamond. 



Fig. 2. (a) Dependence on system size of residues which belong 
to the lowest two two-spinon branches for the Heisenberg model. 
The lowest (second-lowest) branch is indicated by solid symbols 
and 1' (open symbols and 2'). (b) Dependence on TV of the 
poles belonging to the lowest four two-spinon branches (l'-4') 
for a = 0. Symbols have the same meanings with (a). The inset 
shows the dependence on TV of poles at q = n as a function of 
1/TV 2 . 



Fig. 3. (a) Total weight of residues outside the main continuum 
versus q/ir for some values of «(< a c ) and TV. The poles by 
higher-order processes are located outside of the TSC for the 
present values of q (see Fig. 1). (b) Maximum values of the 
same quantity as a function of a for < q/ir < 0.5. The arrows 
near the vertical axes indicate the direction of the dependence 
on TV to the thermodynamic limit around their positions. 



Fig. 4. S(q,ui) for a = 0.2 and TV = 26. The intensity of each 
pole is proportional to the volume of the sphere, and in the same 
scale with Fig. 1. The boundaries of the two-spinon continuum 
for the HS model, &e(q) and Lj u (q), are shown by solid lines. 



Fig. 5. (a) Static susceptibility as a function of q for a = 
0,0.2,0.25,0.4,0.45,0.5 from bottom to top, respectively. Lines 
(except solid line) indicate polynomial fit of 26-site data as a 
guide for the eye. The solid line is the value for the HS model 
given by eq. (3.6). The arrow and cross on the vertical axis indi- 
cate the exact value both for the Heisenberg model and the HS 
model: l/2n 2 . The inset is a further magnification of small-g 
regime to recognize singular behavior of the Heisenberg model, 
(b) Magnified figure of (a) for small q. Symbols are common in 
(a) and (b). 



Fig. 6. (a) K( 2 \q) and (b) K^(q) for some values of a. Symbols 
are common between (a) and (b). Lines (except solid line) are 
guide for the eye. The solid lines in both panels are of the HS 
model given by eqs. (3.7) and (3.8), respectively. 
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Fig. 7. (a) Residues for the MG model. The system with JV = 24 
is shown to include the wave number 7r/2. (b) S(q, u>) for the MG 
model (JV = 26). In both (a) and (b), both contributions from 
\f + and are plotted together, and the weight is indicated by 
the area of the circles. 



Fig. 8. (a) Dependence on JV of residues for the first (1') and 
second (2') lowest dominant branches as a function of q. (b) 
The same for the third (3') and fourth (4') dominant branches. 
Both contributions from and if — are plotted together, (c) 
Dependence on JV of the pole positions corresponding to the 
residues in (a) and (b). The cross and arrow on the right vertical 
axis show the extrapolated pole position of the branch 1' in the 
thermodynamic limit, which means the gap of the MG model: 
A/ J ~ 0.233 in our estimation by the fourth-order polynomial 
fit. The inset shows the movement of the poles at q = tt as 
the system size changes. Since the spectrum is symmetric with 
respect to q = tt/2, we only show 0.5 < q/ir < 1.0. 



Fig. 9. Residues at q = tt as a function of to/ J for some system 
sizes. The connected points with numbers (1' etc.) correspond to 
the poles in Fig. 8. The arrow indicates the peak near the upper 
edge of the main continuum. Since the number of coefficients 
used in the recursion method is different for different systems, 
namely 30 for JV > 22 and 100 for JV < 20, many poles with 
weak intensity appear for the small systems. 



Fig. 10. S(q,u) of JV = 26 for (a) a = 0.35, (b) 0.4 and (c) 
0.45. In this size the ground state is Q = tt. The intensity is 
proportional to the area of the circle, and in the same scale with 
Fig. 7(b). 



Fig. 11. (a) The ground-state energy, (b) the lowest excited en- 
ergy at q = tt/2 and the gap, each as a function of a. Each value 
is extrapolated from the finite-size data (JV = 8-26) by using 
some orders of polynomial fit. The values for the HS model are 
indicated by crosses and arrows on the vertical axes. The energy 
and the gap were previously estimated in some references. |j 



Fig. 12. S(q) for some large values of a. Each line connects the 
data of JV = 26 except for a = 0.5. The arrows along the right 
vertical axis indicate the directions of the dependence on JV of 
S(tt) for JV -» co. 



Fig. 13. Comparison of S(q) between the diagonalization results 
for 5 = and the experimental results for 10K and 20K, which 
are taken from Fig. 2 of ref. 19. The arrows along the right 
vertical axis indicate the directions of the dependence on JV for 
JV -> oo. 



Fig. 14. S(q, ui) at a fixed value of 5 = 0.05 for three values of a, 
(a) 0, (b) 0.2 and (c) 0.4. JV = 26. In each figure, the intensity 
is proportional to the area of the circle, and in the same scale 
with Figs. 7(b) and 10. 



Fig. 15. Lowest excitation at q = tt/2 as a function of 5 for some 
values of a. Shown data are extrapolated from finite-size results 
of JV = 8-24 with using eq. (2.9) and m = 1. The values for 
5 = correspond to Fig. 11(b). The arrow on the vertical axis 
indicates the value of des Cloizeaux-Pearson, tt/2. Shaded lines 
are guide for the eye. The dash-dotted line shows the exactly 
soluble case, as given in Appendix B. 



Fig. 16. (a) Gap of the model eq. (2.1) as a function of 5 for some 
values of a. Each value is extrapolated from the results of JV = 8- 
26 by using eq. (2.9). The arrow on the vertical axis indicates the 
gap for the MG model. If we fit for each a with power functions, 
the resultant powers are 0.71 (a = 0.2), 0.61 (0.35), 0.58 (0.4) 
and 0.51 (0.45). Dashed (dash-dotted) line indicate^ a plausible 
value of a'-NaVoOs for A = 85K and J = 441KEf (A = 98K 
and J = 560KI 49 P ) . The data for small 5 and/or a are omitted, 
because reliable results cannot be obtained by extrapolations. 
Size of the symbols and width of the lines roughly represent the 
inaccuracy of the extrapolation, (b) Magnification of (a) for a 
small-5 regime. Dash-dotted (dashed) line indicates a value of 
0uGeO 3 (A = 24K) for J = 180K (160K). See text. 



Fig. 17. Total energy of the model eq. (2.1) as a function of 5 for 
some values of a. Each value is extrapolated from the results 
of JV = 8-26 by using eq. (2.9). If we fit the data for each a 
with power functions in this range of 5, the resultant powers are 
1.42 (a = 0), 1.30 (0.2), 1.13 (0.35), 1.07 (0.4) and 1.03 (0.45). 
Dash-dotted line indicates the soluble case of 2a + 8 = 1. 



Fig. 18. (a) Dependence on system size of the pole positions of 
the lowest two branches. The parameters roughly correspond to 
CuGeOs. The subtle asymmetry with respect to q = tt/2 is due 
to numerical round-off errors. The inset shows the dependence 
on JV at q = tt for JV = 8-26, fitted by eq. (2.9). (b) Dependence 
on JV of the corresponding residues. The arrow on the right axis 
indicates the extrapolated value for JV — > oo of the branch 1'. 
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Fig. 19. (a) Dependence on N of the pole positions for the low- 
est three dominant branches. The parameters roughly corre- 
spond to a'-NaV205. Concerning the branch 3', we show only 
for q/n ^ 0.4. By '+' we represent poles with tiny weight; to 
avoid mess we plot only for q/n ^ 0.6 and under the branch 3'. 
The inset shows the dependence on N at q = w for N = 8-26 
with polynomial fit. (b) Dependence on N of the corresponding 
residues. 



Fig. 20. Lowest spin excited modes for the ID J-J' model ob- 
tained by the SMA. For the MG case (dash-dotted line), ui = J. 
Otherwise, the diagonalization results are used. Dotted lines are 
guide for the eye. For comparison for a = 0, the exact lower 
(dCP) and upper edges of the TSC are shown with dashed lines. 
The lower edge for a = 0.45 obtained by the recursion method 
(A' = 20-26) is indicated by cross and plus. The SMA and 
the lowest edge for the HS model is indicated by solid and long 
dash-dotted lines, respectively. The extrapolated values of the 
gap (exact) are shown by arrows on the vertical axis for three 
values of a. 



